Kondo effect in coupled quantum dots: a Non-crossing approximation study. 
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The out-of-equilibrium transport properties of a double quantum dot system in the Kondo regime 
are studied theoretically by means of a two- impurity Anderson Hamiltonian with inter-impurity hop- 
ping. The Hamiltonian, formulated in slave-boson language, is solved by means of a generalization of 
the non-crossing approximation (NCA) to the present problem. We provide benchmark calculations 
of the predictions of the NCA for the linear and nonlinear transport properties of coupled quantum 
dots in the Kondo regime. We give a series of predictions that can be observed experimentally 
in linear and nonlinear transport measurements through coupled quantum dots. Importantly, it is 
demonstrated that measurements of the differential conductance Q = dl/dV, for the appropriate 
values of voltages and inter-dot tunneling couplings, can give a direct observation of the coherent 
superposition between the many-body Kondo states of each dot. This coherence can be also de- 
tected in the linear transport through the system: the curve linear conductance vs temperature is 
non-monotonic, with a maximum at a temperature T* characterizing quantum coherence between 
both Kondo states. 
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I. INTRODUCTION 



The recent observations of Kondo physics in 
the electttfflic transport properties of quantum dots 
(QD's)tlEH3cl'El, a tiny_semiconductor box containing a few 
interacting electronsQ, have opened new promising direc- 
tions for experimental and theoretical research of this 
phenomenon, one of the paradigms in condensed matter 
physics. 

The Kondo effect appears in dilute alloys containing 
localized moments as a crossover from weak to strong 
coupling between itinerant electrons of the host non- 
magnetic metal and the unpaired localized electron of 
the magnetic impurity as the temperature iT) is reduced 
well below the Kondo temperature (Tk) .0 Due to spin 
exchange interaction, a many-body spin singlet state is 
formed between the unpaired localized electron and the 
itinerant electrons with energies close to the Fermi energy 
of the metal. This singlet is reflected in the density of 
states (DOS) of the impurity as a narrow peak at low fre- 
quencies: the Abrikosov-Suhl (AS) or Kondo resonance. 
This many-body resonance in the DOS is responsible for 
anomalous properties: In the intermediate temperature 
regime T > Tk, this effect leads to logarithmic correc- 
tions to the magnetic susceptibility x(^)i the linear spe- 
cific heat coefficient 7(T) and the resistivity p(T). Below 
the Kondo temperature, it leads to saturated behavior 
of the magnetic susceptibility x(T) = const, the linear 
specific heat coefficient 7(T) = const and the resistivity 
p{T)-p{0) ~ T2 (Fermi liquid behavior). The Kondo ef- 
fect, being one of the most widely studied phenomenon in 
condensed matter physics, has been studied for decades. 
The first manifestations of the Kondo effect in a linear 
transport property, namely a resistance minimuiB at fi- 
nite temperatures, date back to the early 30 'sE Zero- 
bias anomalies in the nonlinear tunneling conductance. 



the hallmark of Kondo physics, were first observed dur- 
ing the 60's.eI Finally, the first measurements of nonlinear 
transport throuah a single Kondo impurity were reported 
in the mid 90 'sB 

In recent years, spectacular advances in nanotechnology 
have made it possible tcujCjoMa-imentally study Kondo 
physics in quantum dots.BcrBulj These truly impressive 
experiments confirm early theoretical predictions that 
transport through quantum dots in the Coulomb block- 
ade regime shcuild exhibit Kondo physics at low enough 
temperatures.li3 QD's provide the intriguing opportu- 
nity to control and modify the Kondo effect experimen- 
tally: the continuous tuning-|Of the relevant parameters 
governing the Kondo effectQ as well as the possibility 
of studying Kondo physics whe n ih(t , syS it eij i ji .i s . driven 
out of eq.uiljhi'iuYn.. oither by dctiflljilijIiJ't^tSllZl or ac 
voltagesEaOEl'ElaBy, pave the way for the study of 
strongly correlated electron physics in artificial systems. 
Moreover, they provide a unique testing ground in which 
to investigate the interplay of strongly correlated electron 
physics, quantum coherence and non-equilibrium physics. 
More sophisticated configurations of QD's in the Kondo 
regime constitute a growing area of intense investiga- 
tions, both from the theoretica|_aiid_experimental sides. 
Time dependent Kondo physicstjOo'LJ, Kondo physics 
in integer-SBui QD'st3 or QD's embedded in Aharonov- 
Bohm ringg^J are examples of such configurations. 

The study of Kondo physics in mesoscopic or 
nanoscopic systems is not limited to QD's. We can men- 
tion here the recent observation of Kondo phwiics in sin- 
gle atomg22l, moleculesEa, carbon nanotubeg23, scanning 
tunneling microscopy (STM)— experiments of magnetic 
impurities in quantum corralaSil or the anomalous energy 
relaxation in voltage-biased quantma wires and its rela- 
tion to two-channel Kondo physics .t3 

In this paper we will focus on another configuration: 
a system of two coupled quantum dots in the Kondo 
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regime. In view of the recent experimental advances in 
the stu dy , of . q uantum coherence in coupled quantum dots 
(DQD)c3EZlta and the aforementioned studies of Kondo 
physics in quantum dots, it is a timely question to ask 
what happens when a system consisting of two quantum 
dots in the Kondo regime, coupled to each other by means 
of a tunneling barrier, is driven out-of-equilibrium, and 
how the interplay of strongly correlated electron physics, 
quantum coherence and non-equilibrium physics leads to 
new physical scenarios. Previous theoretical studies of 
this problem at equilibrium havs-4«eused on aspects of 

smclcil and on the compe- 



ati-ferromagnetic cou- 



quantum coherence in this syster 
tition between Kondo effect rand, 

pling generated vicuejxchangai^lcaH or via capacitive cou- 
pling between dotscJ. There have hitherto been only few 
attempts to attack this problem in a non-equilibrium sit- 
uation by means of differeat techniques: the equation-of- 
motion technique (EQM)llj, the so-called resonant tun- 
neling approximation^ (valid for T > Tk and equivalent 
to the EOM_method of Ref. |l^ and slave-boson mean 
field theorycj. Here, we present an approach which, for 
the first time, tackles with this non-equilibrium problem 
in a non-perturbative, fully self-consistent and conserv- 
ing way. Our approach is based on a generalizaticp, pf 
the so-called "non-crossing approximation" (NCA)c3'e3 
to the present problem. 

The system of two coupled QD's can be modeled by 
means of two Anderson impurities, each of them coupled 
to a different Fermi sea, and coupled together by means 
of an inter-impurity hopping term. Quantum impurity 
models such as the Kondo and the Anderson impurity 
problems were first introduced in the 60's trying to ex- 
plain the aforementioned anomalous properties of met- 
als in the presence of magnetic impurities. More gen- 
erally, this class of problems and their generalization to 
the lattice constitute one of the paradigms of modern 
condensed matter theory. They typically consist of con- 
duction electrons coupled to sites where there is a strong 
on-site Coulomb interaction, and are believed to describe 
the rich physics of many different strongly correlated 
electron systems like, for example, the heavy fermion 
compounds.Lj Examples of these models are the Hub- 
bard model, the t — J model and the Anderson or Kondo 
lattice models. In all these models the main difhculty 
resides in the fact that usual perturbation theory does 
not apply. On one hand, if the on-site Coulomb repulsion 
exceeds the band width, conventional many-body pertur- 
bation theory in the on-site repulsion does not work. On 
the other hand, the obvious alternative of perturbation 
theory in the kinetic energy is not valid due to the non- 
canonical commutation relations of the field operators in 
the atomic limit. At the heart of the problem is the char- 
acteristic feature of strongly correlated electrons: the dy- 
namics is constrained to a subspace of the total Hilbert 
space. For instance, in the atomic limit of the Hubbard 
model each lattice site can either be empty |0), singly 
occupied | t), I i) or doubly occupied | ti)- The oper- 
ators describing these states and the transitions among 



them, the Hubbard operators, are neither fermions nor 
bosons which precludes the application of usual pertur- 
bation theory (Wick's theorem does not apply). One 
way of circumventing this difficulty is the auxiliary parti- 
cle representation pioneered by Abrikosoit, who first rep- 
resenteid-i local spins b5j-4)seudo-fermionsEij, and later by 
BarnesH and ColemarEj and consists of describing each 
of the states (for each site) as created out of the vacuum 
from the application of a creation operator (bosonic for 
|0) and I ti) and fermionic for | t) and | [) due to quan- 
tum statistics). Each site has to be in one of the four 
states, this is accomplished by constraining the number 
of auxiliary particles to one. Slave particle representa- 
tions allow one to work with usual quantum field theory 
methods provided one works in the constrained subspace 
of the Hilbert space where the number of auxiliary par- 
ticles is one. In particular in the limit of infinite on-site 
interaction, which is case we shall consider in the follow- 
ing, each site can be described by one boson |0) = h^\vac) 



and two fermions | ]) — f}\vac), \ [) = f Avac). This 



particular version of the auxiliary particle representation 
has beenj-termed slave boson (SB) representation after 
Coleman.E3 

Within the SB formulation two non-perturbative ap- 
proaches can be applied to N-fold degenerate Anderson- 
impurity models: 

i) The m|ean-field approximation (MFA) of the slave bo- 
son field j^'Eil only valid for describing spin fluctuations 
in the Kondo regime, correctly generates the low energy 
scale Tk and leads to local Fermi-liquid behavior at zero- 
temperature. The MFA, however, does suffer from two 
drawbacks: a) it leads always to local Fermi liquid be- 
havior, even for multichannel models; b) The MFA has 
a phase transition (originating from the breakdown of 
the local gauge symmetry of the problem) that separates 
the low temperature state from the high temperature 
local moment regime. This later problem may be cor- 
rected by including 1/N fluctuations around the mean- 
field solution.y The generahzation of the SBMFA to the 
present problem, two coupled quantum dots in a non- 
equilibrium situation has been studied in Re£_0 
a) The Non-crossing approximation (NCA)OCJ is the 
lowest order self-consistent, fuUyj-conserving and $ deriv- 
able theory in the Baym sense.Lj It is well known that 
the NCA fails in describing the low-energy Fermi-liquid 
regime. Neglect of vertex corrections prevents from a 
proper description of low-energy excitations. Neverthe- 
less, the NCA has proven to give reliaiile results for tem- 
peratures down to a fraction of Tr-.lZI The NCA gives 
better results in multichannel cases, where the correct 
non-Fermi liauid behavior is obtainedEj. Nonetheless, 
Kroha et aE^ have shown in a series of papers that it is 
possible to develop systematic corrections to the NCA's 
<& functional that cure the low-temperature pathologies 
of the NCA. These systematic corrections (the so-called 
"conserving T-matrix approximation") are able to de- 
scribe Fermi liquid and Non- Fermi liquid regimes on the 
same footing. It is also possible to formulate NCA equa- 
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FIG. 1: Schematic diagram of the double dot system stud- 
ied in this paper. Each dot is coupled independently to one 
lead with couplings Tl and Tu respectively, tc is the interdot 
tunneling term. Note that the role of the inter-dot term is 
twofold: firstly, it generates quantum coherence between the 
two quantum dots; secondly, it establishes a non-equilibrium 
situation, when the chemical potentials are different there is 
a bias voltage across the system hl — fJ^n ~ eV and, then, a 
flow of electrical current through the double dot system. 



tions for finite U by symmetrizing the usual NCA dia- 
grams rwith respect to empty and doubly occupied local 
states.E3 

The generalization of the NCA to time-dependent phe- 
nomena was developed by Langreth and collaborators 
in a series of papersEilH (see also Ref. ^3|). and later 
appUtad—tn ..uijiu-requilibrium transport through jjiuantum 
dotald'EJ'oHEJ and other mesoscopic systems. Ell 
In this work, the NCA is generalized to cope with the 
present problem, namely two Anderson impurities, cou- 
pled to each other by a tunneling barrier, which are in a 
non-equilibrium situation. 

The paper is divided as follows: In Sec. II we formulate 
the Hamiltonian (the general form and its slave-boson 
formulation) which describes the problem. In Sec. Ill we 
briefly review the non-equilibrium Green's function tech- 
nique, real time Dyson equations for the retarded and 
lesser Green's function, that we use in order to formulate 
the problem in its fully non-equilibrium form. In Sec. 
IV we present our generalization of the NCA technique 



to the problem: In Sec. IVA the self-energies obtained 
within our scheme are presented and discussed. In Sec. 
IVB we derive the physical two-particle correlation func- 
tions within the NCA approach. In Sec. IVC we present 
the fermion and boson selfenergies after the projection 
onto the restricted Hilbert space. We present and discuss 
in Sec. V various model calculations for the density of 
states (DOS), both in equilibrium and non-equilibrium 
situations, linear conductance, non-linear current and 
non-linear differential conductance. We give a series of 
predictions for the current and finite voltage differential 
conductance which are relevant for experiments. It is 
demonstrated that measurements of the differential con- 
ductance Q = dl/dV, for the appropriate values of volt- 
ages and inter-dot tunneling couplings, can give a di- 
rect observation of the coherent superposition between 
the many-body Kondo states of each dot. We also give 
predictions for the temperature dependence of the linear 
conductance and for the nonlinear differential conduc- 
tance in the high-voltages regime, where negative differ- 
ential conductance is obtained for low temperatures and 
large inter-dot couplings. An appendix is included to 
discuss the projection procedure used to deal with the 
constraint in the Hilbert space. 



II. MODEL 

A. General formulation 

As we mentioned already, the double quantum dot can 
be modelled as a two-impurity Anderson Hamiltonian 
with an extra term accounting for inter-impurity hop- 
ping. Each impurity is connected to a different Fermi 
sea with chemical potential hl ^ ^ and /i/j — 
respectively {cp = 0)- 



kaeiL.Ry-f ae{L,R},a fcae{t,R},cr cr 

+ ULnL,piL,i + URnR,^nR^l. (1) 

I 



The first two terms in the Hamiltonian represent the 
electrons in the leads and in the dots respectively. In 
these hamiltonians, cj.^^^ ^ (cfei/^.cr) creates (annihi- 
lates) an electron with momentum fc^/i? a-nd spin a in 
the left/right lead, and d^^^j^^ {dh/Ra) creates (anni- 
hilates) an electron with spin a in the left/right dot. 

= efc + Mi/i? = ffc ± ^ and e„a are the ener- 
gies in the leads and the dots, respectively. The third 
term describes the coupling between each dot and its cor- 
responding lead, and determines the coupling strength 



A^'^ = ^^V^Y.k^^^^ ^^'^ ~ neglect the k de- 

pendency of the tunneling matrix element for simplicity). 
The fourth term describes inter-dot tunneling. In the ab- 
sence of inter-dot tunneling, this Hamiltonian describes 
two independent Anderson impurities each of them cou- 
pled to different Fermi seas (typically at different chemi- 
cal potentials). Note that the role of the inter-dot term is 
twofold: firstly, it generates quantum coherence between 
the impurities; secondly, it establishes a non-equilibrium 
situation, when the chemical potentials are different there 
is a bias voltage across the system and, then, there is 
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an electrical current flowing through the double dot sys- 
tem. The last terms describe the on-site electron-electron 
interaction on each dot where ni^i^ „ = j Ra^^L j B.a 
are the number operators for spin a on each dot. The 
on-site interaction parameters are C/^ — I2Cl and 
Un = e^/2Cii where Cl/_r are the dot capacitances. 
The neglect of an interdot electron-electron interaction 
Uinterdotn l-tib) with C/interdot ^ ^^^^ corrcspouds to 
the experimentally accessible limit of small interdot ca- 
pacitance (Clr) as compared with the capacitances of 
each QD to the gates, and implies a vanishiiig| interdot 
antiferromagnetic coupling from this source.E3 Experi- 
mentally, these parameters governing the hamiltonian: 
tunneling couplings, on-site interactions, etc; isan be pur- 
posefully modified by external gate voltagega which al- 



lows to study a variety of rich physical phenomena (spin 
and charge fluctuations regimcj_[H|an|-|equilibrium phenom- 
ena, etc) on the same sample. cIBHcI'EI 

B. Slave-particle representation 

We assume Ul , Ur cx)@, forbidding double occu- 
pancy on each dot. This is a good approximation for 
ultrasmall quantum dots in which the on-site interaction 
is much larger than the coupling strength A^'^ (typi- 
cally more than one order of magnitude). In the limit of 
Ul, Ur oo (i.e, Cl, Cr 0) we can write the Hamil- 
tonian (1) in terms of auxiliary pseudo-fermions and slave 
boson operators plus constraints: 



V 



(2) 



In the slave boson representation, the annihilation oper- 
ator for electrons in the QD's, daa is decomposed into 
the SB operator h\ which creates an empty state and a 
pseudo fermion operator f^a which annihilates the singly 
occupied state with spin a in the dot a: d^cr — > b^afaa 
(d|j^ facrba)- Notc that we have re-scaled the hopping 
parameters Vq — and Vc = -t^- in order to have a well 
defined 1/iV expansion (N being the degeneracy of each 
dot). Finally, the physical constraint is that we must 
work in a subspace of the Hilbert space where the num- 
ber of auxiliary particles (on each dot) is one, namely: 

Qr = = (3) 

As we mentioned before, these two constraints come from 
the physical condition that each dot has to be in one of 
the three states |0), | t) or | |). To simplify the nota- 
tion we consider henceforth that eha — ^Ra = eo- The 
Hamiltonian (2) has two different kind of fermion-boson 
interactions which are given by the vertices in Fig. 2. 



III. GREEN'S FUNCTIONS AND 
SELF-ENERGIES 

At this point, we have reduced the original problem 
described by the Hamiltonian in Eq. (1) to a problem 
of fermions and bosons interacting through the vertices 




FIG. 2: Interaction vertices. Solid, dashed and wavy lines 
represent lead electron, pseudo-fermion and slave boson lines, 
respectively. Each line carries a left (right) index, a) Lead- 
dot hopping vertex V/\/N (Full circle). Tunneling of an elec- 
tron from the left (right) dot to the left (right) lead is rep- 
resented as the decay of the left (right) pseudo-fermion into 
a left (right) slave boson and left (right) lead electron, b) 
dot-dot hopping vertex ^2. (Open circle). Tunneling of an 
electron from the left (right) dot to the right (left) dot is rep- 
resented as the combination of the left (right) pseudo-fermion 
with the right (left) slave boson to decay into a left (right) 
slave boson and right (left) pseudo-fermion. Note that this 
vertex exchanges left and right indexes. 



of Fig. 2 and subject to the constraints in Eq. (3). 
Properties of the physical electrons can be build up from 
the Green's functions of the pseudo-fermions and slave 
bosons (see section IVB). These Green's functions for the 
auxiliary fermions and bosons constitute the basic build- 
ing blocks of the theory. Furthermore, our aim is to study 
the out-of-equilibrium properties of the system; we need. 
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Langreth, see Ref. ^ for a review. 

The time-ordered double-time Green's function are de- 
fined as (sub-indexes are omitted here): 



iG{t,t') EE {TJ{t)p{t')) 
iB{t,t') EE {TMt)b\t')). 



(4) 



' L(R) 



FIG. 3: Diagrammatic representation of the generating func- 
tional $ = $1 + "1>2 of our NCA approximation. Solid, dashed 
and wavy lines represent lead electron, pseudo-fermion and 
slave boson lines, respectively. Each line carries a left (right) 
index. Full circle: Lead-dot hopping vertex V/^/N , Open cir- 
cle: dot-dot hopping vertex ^ . The self-energies are obtained 
by taking the functional derivative of $ with respect to the 
corresponding Green's function, a) Lead-dot functional "l>i 
(leading order 0(1)). b) dot-dot functional $2 (leading order 
0{l/N)). 



then, a fully non-equilibrium description of the dynam- 
ics of the Green's functions of these auxiliary particles. 
The appropriate starting point is to derive equations-of- 
motion (EOM) for the time-ordered double-time Green's 
function of the auxiliary fermion (G) and boson (B) 
fields on a complex contour. A rigorous and well es- 
tablished way to derive ptiiese EOM was first introduced 
by Kadanoff and BaymJiO and has been related to other 
non-equilibrium methods (like the Keldysh method) by 



J 



Here the time ordering operator Tc and the step functions 
operate along a contour c in the complex plane. It will 
not matter in the derivation given here whether c is taken 
to be the Keldysh contour, the Kadanoff-Baym contour, 
or a more general choice. 

The time-ordered Green's functions functions can be 
decomposed in terms of their analytic pieces; 

iG{t,t') = G>{t,t')9{t-t')-G<{t,ty{t' -t) 
iB{t,t') = B>{t,t')e{t~t') + B<{t,t')e{t' -t); (5) 

where G<{t,t') = {fHt')f{t)) and B<{t,t') = {b^t')b{t)) 
are the so-called lesser Green's functions, and G-* {t, t') = 
(fit) Pit')) and B>{t,t') = {b{t)b^{t')) are the greater 
ones. 

The retarded propagators can be written in terms of these 
analytic pieces as: 



iG'^{t,t') = [G>{t,t') + G<{t,t')]0{t-t') 
iB''{t,t') = [B>{t,t')- B<{t,t')]0{t~t') 



The advanced ones can be obtained from G^{t,t') 



(6) 



The basic starting equations follow directly from the Dyson equations in complex time space: 

d 



{t—-eo)Git,t') = S{t-t')+ I dhnt,h)G{h,t'), 



.|5(M') 



S{t-t')+ dtin{t,ti)B{ti,t'). 



(7) 



Applying analytic continuation rules0 we can write Dyson equations in real time space which relate the lesser and 
the greater Green's functions with the retarded and advanced ones: 



(z|-6o)G^(t,0 = 



dh [S]'-(t, h)G< (ti, t') + S^(t, h)G^{ti,t% 
z-B^{t,t') = J dh[U'it,h)B^ih,t') + U^it,h)B''{h,t')] 



(8) 



The retarded (and advanced) Green's functions follow usual Dyson equations 

d 



dt 



eo)G''(i,i') = 5{t-t')+ / dti^^{t,ti)G''{ti,t') 

J — OQ 

i^B^t,t') - 6{t-t')+ / dhW{tM)B^{h.t') 



(9) 



The set of Dyson equations is closed by choosing a suitable approximation for the self-energies E and 11, and hence 
for their analytic pieces S^, S^, 11^ and 11^. We describe in the following section the non-crossing approximation 
used to solve our problem. 



IV. NON-CROSSING APPROXIMATION (NCA) 

A. Self-energies 

We use the NCA techniquell^ for obtaining the 
self-energies {t,t'), l\.Lm){t,t') in Eq. (7) and 



their real time analytic continuations. Hereafter, we 
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FIG. 4: Fermion Self-energy. Solid, dashed and wavy lines 
represent lead electron, pseudo-fermion and slave boson lines, 
respectively. Each line carries a left (right) index. Full cir- 
cle: Lead-dot hopping vertex V/\^, Open circle: dot-dot 
hopping vertex The leading order of this selfenergy is 
0(1/N) +0{l/N'^) 



L(R). 



L(R) 




R(L) 



FIG. 5: Boson Self-energy. Solid, dashed and wavy lines rep- 
resent lead electron, pseudo-fermion and slave boson lines, 
respectively. Each line carries a left (right) index. Full cir- 
cle: Lead-dot hopping vertex V/^/N, Open circle: dot-dot 
hopping vertex The leading order of this selfenergy is 
0(1) + ©(l/iV) 



focus on static non-equilibrium, do voltages, the time- 
translational invariance is thus not broken, i.e all qjiau-. 
titles depend only on the time difference t — i'Jij'O 
Nonetheless, our NCA equations for the self-energies, see 
Eqs. (10-11) below, together with Dyson equations in 
real time Eqs. (8-9) are valid for general situations with 
broken time-translational symmetry by just substituting 
(t — t') {t, t') and solving the fully time-dependent 
problem. 

The generalization of the NCA for time-deoemdent 
phenomena was developed by Langreth et aEdl£3 and 
has been successfully applied, t o non-e quilibrium trans- 
port through quantum dotsEj'r ir ^Mr^, tunnel junctions 
and point contactsEll, noa -tj qy iilibrium dynamics at sur- 
faces and STM studiesBOn. Also, this technique has 
recently been applied to the study of non-equilibrium dy- 



namics in quantum dots in the Kondo regime0 aad to 
the study of non-equilibrium-induced decoherenceEa also 
in quantum dots in the Kondo regime. 

As we already mentioned, this technique can be justi- 
fied as an 1 /N expansion, at lowest order in perturbation 
theory, although it is better regarded as a fully conserv- 
ing, self-consistent, and $ derivable theory in the Baym 
sense.E3 The NCA fails in describing the Fermi liquid 
regime at temperatures much lower than Tk (again, the 
NCA gives reliable results down to a fraction of Tk) due 
to the neglect of vertex corrections in the two-particle 
correlation functions of pseudo-fermions and slave bosons 
(see section IVB). 

Typically, the NCA consists in solving a set of self- 
consistent equations coupling fermion and boson propa- 
gators. To lowest order in both vertices we obtain the 
Baym functional of Fig. 3. This functional consists of 
two terms $ = $i -I- $2- To lowest order in the lead-dot 
vertex we obtain the functional <i>i (Fig. 3. a) which is of 
leading order 0(1) (the order 0{j^) for the vertex is 

not skeleton) . The functional $2 (Fig. 3.b) is constructed 
from the dot-dot vertex and is of leading order 0{j^). In 
principle, it is posible to construct another generating 
functional from the dot-dot vertex which contains off- 
diagonal propagators (this functional can be constructed 
from $2 by replacing all the diagonal fermion and bo- 
son propagators by off-diagonal ones). This functional, 
however, does not contribute to leading order with terms 
0{jf) in the interdot vertex (in other words, there are no 
off-diagonal selfenergies to second order in the interdot 
vertex). The reason being that the bare (in the dot-dot 
coupling sense) off-diagonal propagators are zero. It is, 
thus, consistent to neglect off-diagonal selfenergies within 
the NCA approximation. This way, our NCA guaran- 
tees that all diagrams of leading order 0{^) are included 
within a more general subset of diagrams which includes 
terms to all orders in both vertices. The NCA solution 
obtained from $ = $i+$2 is expressed diagrammatically 
in Figs.^ and |[ These self-energies are obtained by func- 
tional derivation of the Baym functional (Fig. 3) S = ^ 
and n = This guarantees that our approximation 
is conserving. The obtained selfenergies are of leading 
order 0(1) + 0{1/N) (bosons) and 0{1/N) + 0{1/N^) 
(fermions). Applying the real time analytical continu- 



ations of Ref. 71 to the complex-contour-time-ordered 
fermion self-energy in Fig. ^ we obtain the lesser, greater 
and retarded components: 
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■N 



At-t')B>{i'-t)Bt,m{t-t')} 
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RiL), 



(10) 



Where KL(^ji)^cr{t — t') is the fermion propagator in the left (right) lead. The corresponding expressions for slave- 
boson self-energies (Fig. 5) are: 



n< 



n> 



n 



HR) 



it-t') 



N 
1 

N 



^ Y.i^LiR)Jt' - t)Gl(R).At - t') + Kl(R.)Ai' - ^)G<(^),.(i - t')} 

(T 



BR{L){t - t')GR(^L)A'^' ^ ^)Gt{R)A'^ - t') 
BR(L)i* - t')GRi^L)A^' ^ ^)Gt(R)^ai^ - t')}- 



(11) 



Eqs. (10) and (11) are the unprojected full NCA self-energies coming from the generating functional. The projection 
of these quantities onto the physical subspace Qa£{L,R} = 1 is discussed in Appendix A. 



L(R) 




L(R) 

FIG. 6: Diagrammatic representation of the physical two- 
particle correlation function within the NCA-approximation. 
The neglected vertex corrections are 0{^)^ 



B. Physical correlation functions 

The physical lesser and greater correlation functions 
[a e {L, R}) are: 



A<,{t-t') S (dL(t')da.(i)) 



In terms of slave operators they become the two-particle 
correlation functions: 

A<,it-t') EE {fUt')bc.{t')hi{t)U{t)) 
A>^{t-t') EE {bUt)U^{t)fUt')b^{t'))- (13) 

The evaluation of these two-particle correlation functions 
would require in principle a further diagrammatic expan- 
sion. Within the NCA approximation, however, one ne- 
glects vertex corrections^ and keeps only the lowest or- 
der term in the expansion of the two-particle correlation 
function (Fig. 6): 



A^t ' t') 



{fiAt')fUt)){ho.{t')bl{t)) 
G<At~t')B>{t'-t) 

{bUt)b^{t')){UAt)fL{t')) 

B<{t'-t)G>Jt-t'). 



(14) 



(12) 



Using the identities 



G 



LiR)jt-t') = ^[GiiR)At-n 



R> 

^L{R),a\ 



L(R)A^ ^ ^')] - G'l^R^ „{t - t') 



Eq. (14) can be rewritten as: 

A<,{t-t') = G<At~t'){t[B:,At' 
A>,it-t') = B<it'^t){z[GlAt- 



Gl,...,^.. 
Bl(r)A^ - i')] + B^R-^^t 



t') 



-t)-Bl,{t' -t)]+B<,{t' -t)] 
t')-Gl,At-t')]-G<,,{t^t')]. 



(15) 



(16) 
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Now, according to the projection procedure explained in appendix A (see also Ref. ^ the terms G'^^{t—t')B^^{t' — t) 
and B^{t' — t)G^ „{t — t') have to be projected out (they are of order C'(e~^''^'*'° )). This is accomplished by making: 

Bl(Huii-i') - ABUR)At-^')-BliR)At-t')] (17) 
which gives the following physical correlation functions: 

AUt-t') = ^B<{t'-t)[G:Jt-t')-GlJt~t')] 

A:j:\t t') ^ GZ\t - t')B<{t' -t)- G<At - t')BZ\t' - t\ (18) 

These Green's functions have to be calculated with the corresponding projected selfenergies, as discussed in the next 
section. 

C. Physical selfenergies 

The final set of projected selfenergies is (see Appendix A for details): 



S 



,(r) = (^j^,>(H)^.(r)+z(|)^[G^(^),,(r)-G°(^)^,(r)] ^y'^ x^£(^)(r) 



R(L) 



^UR)Ar) ^ |^i^^(^),.(r) + »(^)' ^Xw"^ [<L)(-r) ^ <L)(-r)]| X g<(^,)(r) 
n^wW = (^E^L%),.(--) + *(^)'E^|^['5«(L),.(--)-G^^^^^ (20) 



GR(L),ai^'^^ 



Where we have introduced the notation t — t — t'. The 
conduction electron propagators K are defined in terms 
of the Fourier transforms of the bare conduction elec- 
tron propagators (namely, without dot-lead coupling) as: 

Kl^ie) = 2^Y.k^VH{e~ ek^)f^{e), where /<(e) = 
e,j(.-L)+i is the Fermi function and /> (e) = 1 - /< (e) 
(see Ref. |6^). The Green's functions G, B do not include 
inter-dot hopping. Finally, the factors and Zji can be 
identified with the left and right charges in the absence 
of inter-dot hopping. They can be obtained from the 
left and right charges of two independent single impurity 
problems (at different chemical potentials ^.l and /i^ re- 
spectively) . It is important to emphasize two aspects of 
the projection: i) the simplification of the propagators 
[K K , G G and B B) is required by the projec- 
tion procedure (see Appendix A) and is not an additional 
approximation; ii) this should not be construed to imply 
that there is no inter-dot correction in the slave-particle 
Green's functions that enter into the physical correlation 
functions of Eq. (18). 



Eqs. (19-20) constitute the main result of this section. 
The projected self-energies inserted in the appropriate 
Dyson equation, give an overall result in Eq. (18) for 
the physical correlation functions which has the correct 
order. Of course, in the absence of inter-dot coupling we 
recover from Eqs. (19-20) two independent sets of NCA 
equations for the left and right single impurity problems. 
These equations are in agreement with the ones previ- 
ously obtained in Refs. [l^ , |6l] , |67| . 



The equations for the self-energies, together with the 
Dyson equations for the retarded and lesser propagators 
and the normalization conditions, close the set of equa- 
tions to be solved. We numerically iterate them to con- 
vergence. 
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FIG. 7: Equilibrium density of states (DOS) for different val- 
ues of the inter-dot hopping tc ~ 0.0, 1.0, 1.2, 1.4, 1.6, 1.8, and 
2.0. The curves are shifted vertically for clarity, a) Full DOS. 
The splitting in the DOS corresponds to the formation of 
bonding and anti-bonding combinations of the single parti- 
cle levels due to inter-dot tunneling, b) Blow up of the low 
frequency region around the Kondo peak. As the inter-dot 
coupling increases, the Kondo peak also splits. Importantly, 
this splitting, which is a manifestation of quantum coherence 
between the two many body Kondo states on each dot, is 
much smaller than the splitting of the broad peak, see main 
text. 



V. RESULTS 
A. Density of states 

Here we present results for the left and right dot den- 
sities of states (DOS), both for equihbrium and finite 
voltage {fiL = V/2, = —V/2) situations. We 
use the following parameters in the calculations (un- 
less otherwise stated): eo = —2.5, T=0.003, all energies 
are given in units of r(e) = 2A(e) = '^'^VoJlk^i'^ ~ 
Efc) — ~ ^fc)- Each lead is described by a 

parabolic density of states centered at the chemical po- 
tential and with a bandwidth W = 2D ~ 12. The 
Kondo temperature corresponding to these parameters is 
r]^ ^ 3.7 X 10~^ (here, the superscript "0" means with- 




u 0.2 
0.1 








-40 -20 20 40 
(B( Units of T°j,) 



-40 -20 20 40 
Q)( Units of T°j^) 

FIG. 8: Temperature dependence of the equilibrium density of 
states (DOS) around the Fermi level. The arrows indicate the 
direction of increasing temperatures. T% = 0.00037. Main 
figure: DOS of the coupled dot problem for tc = 1.6. The 
structure originating from the inter-dot coupling is still visible 
at temperatures T > lOT^. The inset shows the reduction 
of the Kondo peak for the single impurity problem as the 
temperature increases. At the highest studied temperature 
T ~ 20T^ the Kondo peak is almost suppressed as compared 
with the coupled dot system. 



out inter-dot coupling, namely the Kondo temperature of 
the single impurity problem corresponding to these pa- 
rameters), as calculated from the Bethe ansatz analytical 
solution {N = 2)mi: 



rpO 



r(l + l/2)D^(2A/7rD^)i/2e-'^l'^ol/2A^ (31) 



r(.x) is the gamma function and the re-scaling Dj. — 
g-i/2£) accounts for the assumed parabolic DOS in the 
leads instead 0|Lthe rectangular one used in the Bethe 
ansatz solutionEll. Note, finally, that in order to com- 
pare with the slave-boson mean field (SBMF) results 

It is known from slave-boson MFAeilcSO and from nu- 
merical renormalization groupEJ calculations that the ef- 
fective Kondo temperature T^'^ of the double dot system 
grows exponentially with the inter-dot hopping. Unfor- 
tunately, it is not possible to extract an analytical expres- 
sion for the Kondo temperature from our set of coupled 
NCA equations. We choose then relatively high temper- 
atures T > T]^ in all our calculations in order to prevent 
the expected low-temperature pathologies should the ef- 
fective Kondo temperature increase with the inter- 
dot hopping. 

In Fig. 7, the QD density of states (DOS) at equilibrium 
(here, of course, the left and right dot are equivalent) is 
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FIG. 9: Non-linear transport properties of the DQD system 
for different inter-dot couplings for T = 0.003 ^ ST]^. a) 
Current-voltage characteristics, b) Differential conductance 
at finite voltage. The zero-bias anomaly first broadens and 
then splits with increasing inter-dot hopping. The splitting 
of the zero bias anomaly reflects quantum coherence between 
the two many-body Kondo states on each QD. 



plotted for increasing values of the inter-dot tunneling. 
The full DOS (Fig. 7a) shows the splitting of the main 
peak (energy scale for charge fluctuations) originating 
from the inter-dot coupling which generates quantum co- 
herence between the dots. The splitting in the DOS cor- 
responds to the formation of bonding and anti-bonding 
combinations of the single particle levels, i.e. e± = eoT^c 
due to inter-dot tunneling. Fig. 7b, shows a blow up of 
the low frequency part of the DOS around the Fermi 
level. As we increase the inter-dot coupling, the Kondo 
peak also splits into bonding and anti-bonding combina- 
tions. Importantly, the energy scale for this splitting of 
the Kondo peak, which is a manifestation of quantum 
coherence between the two many body Kondo states on 
each dot, is much smaller than the one corresponding to 
the splitting of the broad peak (which is a manifestation 
of coherence between single particle states). We have, 
then, Ae = 2ic << Ae = 2tc, where Ae and Ae = 2tc 
are the splitting of the Kondo peak and the single parti- 
cle splitting, respectively. This reduction of the splitting, 
namely tc << tc, is caused by the strong Coulomb re- 
pulsion on each dot. Typical values of this splitting are 
in the range ~ lOT^ — 40T^ (note that the single parti- 
cle splittings are in the range ~ lO^T^ - lO'^T^). These 
obtained values for the reduced splitting of the Kondo 
resonance are in good semiquantitati ve . agrocni Lcnt with 
the mean-field slave boson calculation.c2fEJc3o 



The behavior at different temperatures is studied in 
Fig. 8 where we plot the DOS of the coupled dot prob- 
lem (Fig. 8a) for tc = 1.6 at different temperatures. 
The splitting originating from the inter-dot coupling is 
still visible at temperatures T > lOT^. For comparison, 
we show in Fig. 8b the reduction of the Kondo peak for 
the single impurity problem as the temperature increases. 
At the highest studied temperature T ~ 20T^ the Kondo 
peak for the single impurity system is almost suppressed 
as compared with the coupled dot system. This is in good 
qualitative, agr ti om ftnt with the previous statement that 
T^^ > Tj^Eiyy-y. it is worth noting that the sphtting 
of the Kondo resonance is robust at temperatures higher 
than T^; experimentally this is of the most relevance: 
according to this result, the experimental conditions for 
studying Kondo physics in coupled QD's are less demand- 
ing than in single QD's (temperatures much lower than 
are needed in order to observe Kondo- trdatcd-features 
in the transport properties of single QD'soaoLttl). 



B. Non- linear transport properties 

We have proven in the previous section that the inter- 
dot coupling generates quantum coherence between the 
dots. This quantum coherence is reflected in the DOS 
of each QD as a splitting, both in the charge fluctuation 
and spin fluctuation parts of the spectrum. We are inter- 
ested in Kondo physics and the obvious question we want 
to answer is thus: Can we observe the splitting of the 
Kondo peak, induced by the inter-dot coupling, in a dif- 
ferential conductance measurement? The answer to the 
previous question is non-trivial because we are dealing 
with the non-equilibrium physics of strongly correlated 
electrons and hence the spectral functions are expected 
to strongly depend on the applied bias voltage (shift and 
broadening of the peaks). In other words, the differential 
conductance curve does not just mimic the zero-voltage 
DOS (as it does for non- interacting electrons). From the 
experimental point of view this is a timely and crucial 
question: the observation of such a splitting would prove 
the remarkable phenomenon of quantum-coherence be- 
tween the two Kondo many-bodM-.states on each dot. Ex- 
periments by Oosterkamp et ac3 and Blick et oH have 
proven quantum coherence between single particle states 
in coupled QD's. Also, some signatures of coherence be- 
tween Kondo states in a double quantum dot system have 
been reported recently by Jeong et al in Ref. |3^. The 
first step in order to answer our question is to calculate 
the current through the double dot system. We follow the 
standard non-equilibrium approach to transport through 
a region of interacting electronsliilEj and relate the cur- 
rent through each dot to its retarded and lesser Green's 
functions: 

/ae{L,fl}--^ J Ar(e)[2W^(e)/4e) + A<(e)]. 

(22) 
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FIG. 10: Non-equilibrium DOS {^l = V/2 = 12.57^, 
/Ltfl = —V/2 = — 12.5Tx) for different values of the inter-dot 
hopping tc = 0.0, 1.0, 1.2, 1.4, and 1.6. The curves are shifted 
vertically for clarity. Top: left DOS. Bottom: Right DOS. 
The arrows mark the position of the chemical potentials. 
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FIG. 11: Non-equilibrium DOS at tc = 1.6 for different volt- 
ages in the range V = lOTg- through 50T^. Top: left DOS. 
Bottom: Right DOS. The arrows mark the directions of in- 
creasing voltages. 



Here, j4^(e) and A^{e) are the Fourier transforms of 
the retarded and lesser physical Green's functions of Eq. 
(18). The total current through the system is calculated 
as 7 = i^^^^j^. The differential conductance Q = dl/dV 
is calculated by numerical differentiation of the current- 
voltage (I-V) curves. 

We study in Fig. 9 the non-linear transport proper- 
tics of the DQD system. Wc plot in Fig. 9a the I-V 
characteristics for different values of the inter-dot hop- 
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FIG. 12: Non-linear transport properties of the DQD system 
{tc ~ 1.6) for different temperatures as a function of the 
applied bias voltage, a) Current-voltage characteristics, b) 
Differential conductance at finite voltage. 



ping. As the inter-dot hopping increases, the low-voltage 
differential conductance grows while at the same time 
starts to deviate from an Ohmic behavior. At large volt- 
ages the current saturates, the differential conductance 
nears zero and even becomes slightly ncigativc for the 
largest tc- These features are better brought out in a 
plot of the differential conductance at finite voltage (Fig. 
9b). As we increase the inter-dot hopping, the zero bias 
Kondo anomaly broadens and splits. We can attribute 
this broadening to the aforementioned increase of the 
effective Kondo temperature as a function of the inter- 
dot hopping. For large inter-dot tunneling couplings the 
zero-bias anomaly splits. The splitting of the zero bias 
anomaly is an unambiguous indication of quantum co- 
herence between the Kondo states on each dot. 

In Fig. 10 we plot the non-equilibrium DOS {hl = 

V/2 = 12. 5T]^, HR = -V/2 = -12.^T%) for the left (top 
figure) and right (bottom figure) coupled quantum dots. 
For the uncoupled situation [tc = 0) each DOS has a 
Kondo peak around each chemical potential as expected. 
With increasing inter-dot hopping, the behavior of each 
DOS becomes quite complex. The Kondo peak on each 
side splits into two peaks while at the same time the 
whole spectral weight near the Fermi level shifts to lower 
frequencies. Furthermore, these split peaks are asym- 
metric, they have different heights and spectral weights 
(it is important to mention here, however, that the NCA 
is known to overestimate the asymmetry of the peaks 
because it does incorrectly treat potential and spin flip 
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FIG. 13: Comparison of the temperature dependence of the 
hnear conductance in a single quantum dot (solid line) and a 
double quantum dot with tc = 1.6 (dashed line). The linear 
conductance for the single dot follows the usual logarithmic 
increase at intermediate temperatures followed by a satura- 
tion near the unitary limit. The linear conductance for the 
double dot case shows a nonmonotonic temperature depen- 
dence, it increases for decreasing temperatures in the region 
T > T* whereas it decreases in the region T < T* . The tem- 
perature scale T* characterizes quantum coherence between 
both dots in the Kondo regime. Note that in order to compare 
with the single dot case the temperature has been scaled with 
respect to the Kondo temperature of the single dot problem, 
see main text. 



scattering on equal footin^Sl). As the intcr-dot hopping 
increases, the lower (upper) band of the left (right) DOS 
moves to lower (higher) frequencies, while increasing its 
height, until it matches with the upper (lower) band of 
the right (left) DOS. As an example, for tc = 1.6 the 
lower peak on the left DOS and the upper peak on the 
right DOS approximately match at w ~ /i^. As a re- 
sult, there is a peak in the differential conductance at 
V — 25T^ for this inter-dot couphng. 
Also interesting is to study how the DOS evolves as a 
function of the applied voltage for a fixed inter-dot cou- 
pling. This non-trivial behavior of the DOS versus ap- 
phed voltage is studied in Fig. 11 where we plot the 
non-equilibrium DOS for tc — l-G and different voltages 
from V = lOT^ to V = 50T^ in intervals of AV = 5T^. 
As the voltage increases, the left (right) DOS moves to 
higher (lower) frequencies such that the middle point be- 
tween the split Kondo peaks lies approximately at the left 
(right) chemical potential (this discussion is only quali- 
tative, note that even for the uncoupled case the Kondo 



peaks do not lie exactly on each chemical potential). 

The temperature dependence of the current and dif- 
ferential conductance are plotted in Fig. 12. Several 
features in these curves are noteworthy. If we focus first 
in the differential conductance (Fig. 12. b) we see that 
the splitting of the zero-bias anomaly can be resolved for 
temperatures T < lOT^. For higher temperatures the 
splitting can no longer be resolved and, instead, a broad 
zero-bias anomaly is obtained. Also important to men- 
tion is the non-monotonic behavior of the linear conduc- 
tance Q = dl /dV\v=o with temperature. Starting from 
high temperatures, the linear conductance first increases 
for decreasing temperatures, indicating the appearance of 
Kondo physics. This behavior saturates at the tempera- 
ture for which the splitting is resolved (here T ^ lOT]^) 
and then the linear conductance decreases for decreas- 
ing temperatures. This behavior can be easily explained 
by noting that the linear conductance at finite tempera- 
tures is a convolution of the DOS around the Fermi level 
with the derivative of the Fermi function (whose full- 
width at half-maximum is 3.5T). When the width of the 
derivative of the Fermi function is smaller than the split- 
ting of the Kondo peak this convolution is very small, 
due to the small spectral weight around the Fermi level 
when the Kondo peak splits, explaining why the linear 
conductance decreases when lowering the temperature. 
This non-monotonic temperature behavior is an indirect 
proof of the formation of the splitting (in single dots in 
the Kondo regime the linear conductance monotonicaly 
increases, until it saturates in the Fermi liquid regime, 
for decreasing temperatures). We show this behavior in 
Fig. 13 where we compare the temperature dependence 
of the linear conductance of a single quantum dot (solid 
line) with the temperature dependence of the linear con- 
ductance of a double quantum dot with tc — 1.6 (dashed 
line). The linear conductance for the single dot follows 
the usual logarithmic increase at intermediate temper- 
atures followed by a saturation near the unitary limit. 
The linear conductance for the double dot case shows a 
nonmonotonic temperature dependence, it increases for 
decreasing temperatures in the region T > T* whereas it 
decreases in the region T < T* . The temperature scale 
T* (which is the temperature for which the splitting is 
resolved in Fig. 12. b) characterizes quantum coherence 
between both dots in the Kondo regime. Note that in or- 
der to compare with the single dot case the temperature 
has been scaled with respect to the Kondo te mperature 
of the single dot problem, Tk = D ^2 A/ ''^''"^/^^ . 
Finally, it is important to mention here is that the NCA is 
known to overestimate the Kondo peak amplitude (and 
then the linear conductance) when calculated from the 
density of states^ Typical overestimates are within the 
range 10 — 159cI1j. Keeping this overestimation in mind 
(which for temperatures T < 4 x IO^^Tk leads to an 
overshooting of the unitary limit in the single dot case. 
Fig. 13 solid line), we purposefully show results at low 
temperatures where the temperature dependences of the 
linear conductance for single and double dot cases com- 
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FIG. 14: Non-linear transport properties of the DQD sys- 
tem {tc = 1.6) for two different temperatures as a function 
of the apphed bias voltage, a) Current-voltage characteris- 
tics, b) Differential conductance at finite voltage. At large 
voltages the system develops regions of negative differential 
conductance accompanied by fluctuations in the current. We 
speculate that these fluctuations could originate from a dy- 
namical instability in this region of voltages, c) Blow-up of the 
low voltage region in the differential conductance. The extra 
structure at low voltages (small zero-bias anomaly-|-satellites) 
is originating from the splitting due to the applied voltage. 



pare best. 

Finally, we comment on the temperature dependence 
of the differential conductance at large voltages (see Fig. 
12. a). At low temperatures the slope of the I-V charac- 
teristics at large voltages approaches zero and eventually 
becomes slightly negative, namely the I-V characteristics 
present negative differential conductance (NDC), at the 
lowest temperatures. The slope of the differential con- 
ductance increases gradually as one increases the temper- 
ature. For the highest temperature studied (T — 20T^, 
dashed-dotted line) no traces of NDC are found even for 
very large voltages. 

We study this NDC behavior in Fig. 14, where we com- 
pare the I-V characteristics (Fig. 14a) and differential 
conductance (Fig. 14b) of the system at high {T = 20T^) 
and low temperatures {T = 5T]^). For the low temper- 
ature situation, the slope at large voltages does indeed 
develop NDC for voltages V > 50T^. At V ~ 50T^ the 
differential conductance becomes zero and the current 
smoothly starts to decrease as one increases the dc volt- 
age. However, the situation changes drastically at larger 
dc voltages where our numerical results for the current 
rapidly develop a wiggly pattern. The appearance of this 



fluctuating pattern in the numerics is accompanied by a 
breakdown of current conservation, namely Jl = —Jr is 
no longer fulfilled. We do not have a conclusive answer 
for the appearance of this pattern in the current, but, 
nonetheless, the fact that it appears in the NDC region 
together with a breakdown of the condition Jl + Jr. = 0, 
allows us to speculate that it may be reflecting the for- 
mation of a dynamical instability where a time dependent 
current spontaneously develops in response to the static 
applied voltage (with a non-zero displacement current 
Jd^sp{t) = -(Jl + Jr)/2 = -f wWe Qdd is the 

charge accumulated in the double dot). 113 Although the 
analysis of time-dependent phenomena is beyond the pur- 
pose of this work (our numerical scheme is only valid for 
time-translational invariant situations), we just mention 
that this kind of dynamical instability, rather typical in 
non-linear systems presenting NDCO, have been-cecently 
reported in single QD's in the Kondo regime .lJlSEZI It is, 
thus, reasonable to expect similar dynamical instabilities 
in DQD's, which motivates our speculation. 

We finish this part with two remarks. The first is that 
this NDGJias been previously reported in the context of 
SBMFT.E] Importantly, the NDC features obtained here 
are smooth (the dl/dV evolves from zero to negative val- 
ues in a smooth manner) and gradually dissapear as the 
temperature increases (as already anticipated by us in 
Ref. ^) in contrast with the ones obtained within the 
SBMFT (sharp transitions between the high and low cur- 
rent regions). These sharp transitions can be attributed 
to the lack of fluctuations (ouantum and thermal) of the 
boson fields in the SBMFT.El 

The second remark is that the low-voltage part of 
the differential conductance curve at the lowest temper- 
ature (Fig. 14c) does also develop new fine structure 
(extra peaks). The differential conductance develops a 
small zero-bias anomaly and satellites separated from 
zero at AV ~ ±10T^. These new structures in the 
differential conductance are in agreement with, the , ones 
previously reported in two-level quantum dotsriri and 
coupled ouantum dots in the limit of strong inter-dot 
repulsiono and can be attributed to the extra splitting 
induced by the applied voltage: the voltage splits the 
peaks in the left and right spectral functions, a peak 
in the differential conductance occurs when these split 
peaks cross each other. The agreement in only quali- 
tative though. In Refs. such crossings occur at 
Ae — V where Ae (a fixed quantity) is the energy separa- 
tion between single particle levels in the two-level quan- 
tum dotCj'Ej or the the energy separation between the 
bonding and-anti-bonding levels in the coupled quantum 
dot system.Ej On the contrary, the peaks in the differ- 
ential conductance of our calculation appear at much 
lower frequency scales. As mentioned before, our cal- 
culation includes the strong renormalization of the lev- 
els due to electronic correlations and due to the volt- 
age. The crossings, hence, appear at voltages for which 
Ai{V) = V (namely, e+ + V/2 = e_ - V/2), where 
Ae(y) = e_ — e_|_ = 2tc{V) is the voltage-dependent en- 
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FIG. 15: Non-equilibrium full DOS at low temperature 
(T = 5T^) and tc = 1.6 for different voltages V = 
25T^, 50T^, lOOT]^. The applied voltage induces extra split- 
tings in the bonding and anti-bonding combination of the 
Kondo peak. As a result, four peaks can be clearly resolved 
in the full DOS at high voltages (these peaks are marked with 
asterisks for the highest voltage in the figure). 



ergy separation between the anti-bonding and bonding 
combinations of the Kondo peak (which, again, is much 
smaller than the single-particle splitting 2tc)- 

Fig. 15, where we plot the full spectral function at a 
finite voltage, illustrates this phenomenon. Each peak 
splits by zkV/2. As a result the full DOS develops four 
peaks, the combinations and e_ ztV/2, that can 

be clearly resolved at high enough voltages. These split 
peaks are marked with asterisks for the highest voltage 
in the figure, the distance between consecutive peaks is 
twice the renormalized inter-dot hopping, the distance 
between alternate peaks is the voltage. We mention, in 
passing, that the observation of this fine structure in the 
differential conductance would constitute a direct proof 
of the voltage-induced splitting of the Kondo resonance. 
Here, the splitting associated with the inter-dot hopping 
serves as a testing tool, similarly to that provided by 
an external magnetic field in single quantum dotsE3 (the 
quantity Ae of our previous discussion being now the 
zeeman splitting in a single quantum dot with an exter- 
nal magnetic field) to check the voltage- induced splitting. 
Different proposals for measuring this voltage- induced, 
splitting are the subject of current active researchE30'E3 

We support our previous paragraph by studying the 
temperature dependence of the non-equilibrium full DOS 
at F = IOOTO. (Fig. 16). At high temperatures, T > 2ic, 
the splitting coming from inter-dot coupling can not be 
resolved and the coupled dot system is equivalent to a sin- 
gle dot with a broad Kondo peak (coming from a convo- 
lution of the bonding and anti-bonding peaks with ther- 
mal broadening) . The width of this effective Kondo peak 
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FIG. 16: Temperature dependence of the non-equilibrium full 
DOS (V^ = 100r£^) and tc = 1.6. At high temperatures, the 
splitting coming from inter-dot coupling can not be resolved 
and the coupled dot system is equivalent to a single dot with 
a broad Kondo peak coming from a convolution of the bond- 
ing and anti-bonding peaks with thermal broadening. This 
eff'ective single Kondo peak is split by the voltage as expected 
(vertical marks). Further lowering of the temperature allows 
the resolution of the inter-dot-induced splitting (asterisks) . 

is thus larger than 2tc- As expected, a finite voltage, 
V > T, splits this effective single Kondo peak into two 
peaks separated by V (Fig. 16, thick solid line). Fur- 
ther lowering of the temperature allows to resolve the 
inter-dot-induced splitting resulting in extra peaks su- 
perimposed to the ones induced by the voltage (Fig. 16, 
asterisks). 

We finish by commenting on the observability of the ef- 
fects described in this section: We obtain in our calcu- 
lations splittings in the differential conductance of the 
order of ~ 50T^. Typical Kondo temperatures in quan- 
tum dots are of the order of a few fj.eV (for instance, the 
Kondo temperature is ~ 4 — in the experiment 

of Rcf. which gives splittings well within the resolu- 
tion limits of state-of-the-art techniques (remember that 
lfj.eV - lOmK). 



VI. CONCLUSIONS 

We have theoretically studied the transport proper- 
ties, both equilibrium and out-of-equilibrium properties, 
of a coupled quantum dot system in the Kondo regime. 
We have modeled the double quantum dot system by 
means of a two-impurity Anderson Hamiltonian with 
inter-impurity hopping and infinite on-site interaction on 
each dot. The Hamiltonian, formulated in slave-boson 
language, is solved by means of a generalization of the 
non-crossing approximation (NCA) for the present prob- 
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lem: two quantum dots in the Kondo regime, coupled to 
each other by a tunnehng barrier and with an apphed 
voltage across them. We have provided benchmark cal- 
culations of the predictions of the non-crossing approxi- 
mation for the linear and nonlinear transport properties 
of coupled quantum dots in the Kondo regime. We give a 
series of predictions that can be observed experimentally 
in linear and nonlinear transport measurements through 
coupled quantum dots in the Kondo regime: 

i) The nonlinear differential conductance G — dl/dV 
directly measures the transition (as tc increases) from 
two isolated Kondo impurities to a coherent bonding 
and anti-bonding superposition of the many-body Kondo 
states of each dot. For increasing inter-dot couplings, the 
zero-bias anomaly first broadens and then splits. The 
later case corresponds to transport which is optimized 
for a finite bias voltage matching the splitting between 
these two bonding and anti-bonding states. 

ii) The effective Kondo temperature of the coupled sys- 
tem increases with the inter-dot coupling. This is re- 
flected as broadening of the zero-bias anomaly. 

iiij The non-monotonic temperature behavior of the lin- 
ear conductance Q = dl /dV\v=o is an indirect proof of 
the formation of the splitting. Starting from high tem- 
peratures, the linear conductance first grows for decreas- 
ing temperatures, indicating the appearance of Kondo 
physics. This behavior saturates at the temperature for 



which the splitting is resolved. Further lowering of the 
temperature produces a decrease of the linear conduc- 
tance: The curve linear conductance vs temperature has 
a maximum at a temperature scale T* characterizing 
quantum coherence between the two quantum dots. 

iv) The differential conductance at large voltages can be- 
come negative (NDC). We speculate that the system can 
develop dynamical instabilities around this region. 

v) At low enough temperatures, it is possible to resolve 
extra structures in the differential conductance coming 
from the splitting induced by the applied bias voltage. 
We hope our work will inspire and encourage experimen- 
tal investigations of Kondo physics in coupled quantum 
dots and related systems. 
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I 

APPENDIX A: PROJECTION 



Here we discuss the evaluation of operator averages within the restricted subspace of the Hilbert space with the 
constraints of Eq. (3). The formal expression for the expectation value of an operator in this subspace can be written 
as: 

(6)q,=i,q«=i = ^ Tr{e-P^"o-^^-^.~^^-Nn)SQ^^,SQ^,^Tc[Sci~^, oo)6]}, (Al) 

^Ql = 1,Qr = 1 

where Tc orders operators along a complex contour, the hopping terms are treated as perturbations (i.e Hq contains 
the isolated regions of the problem, leads and dots, before they are connected), and the partition function is given by: 

^Q,=i,Q«=i = rr{e-'3(««-^-^--'^«^«)<5Q,,i<5Q,,iTc[^c(-oo, oo)]}. (A2) 

Since the charge operators commute with the Hamiltonian each constraint can be incorporated by a Kronecker delta 
bmction in the statistical averages of Eqs. (A1-A2). To relate averages in the restricted ensemble with the ones 
corresponding to an. imrestricted ensemble we represent each Kronecker delta function as an integral over a complex 
chemical potentiatJo (see also appendix D in Ref. ^ and chapter 7 of Ref. 0): 

<5q«,i = ^ rdXne'^^^-(Q-''\ (A3) 



{0)Ql = 1,Qr = 1 - 



,=1 '27r 



' dXL r dXRe'^^^e'^^^ZGc{d}GC- 



(A4) 
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This way, we can relate the average in the constrained ensemble with an average in the grand canonical ensemble 
which can be written as: 

Zgc 

Zgc = Tr{e~'^^"'>-^''■^'^-^"'^''+'^'■^'■+'^''^''^Tc[Sc{-oo,oo)]}. (A5) 

This average inside the integral in Eq. (A4) now obeys a linked cluster theorem and we can use conventional field 
theory. In principle, we can stop here, evaluate the averages in the grand canonical ensemble and projecting to the 
physical ensemble by a final integration over the chemical potentials. Further simplification can be gained, however, by 
noting that the grand canonical partition function Zqc can be rewritten as a sum over canonical partition functions: 

Zgc^Y. E ^c(QL,Qfl)e-^^'^'^^e-'^^^-««; (A6) 

Qi=OQn=0 

and by expanding the expressions in the grand canonical ensemble as power series: 

Zgc = ^c(0, 0) + Zc{l,G)e~'^^^ + Zc{0. l)e-'^^« + Zc(l, l)e~-P^^ e"^^'' + ... 
{0)gc = (0)°^° + (0)''°e-*''^^ + (d)"'ie-*''^« + (C')'''e-'^^^e-'^^« + ... (A7) 

Inserting these power series expansions inside the integral in Eq. (A4), the only terms that survive are: 

(0)q,=i,q«=i = — ^[Zc(0,0)(6)i'i + Zc(l,0)(d)"'i + Zc(0,l)(6)i'°], (A8) 
^cU, ij 



where we have used (O)"^'^ = which is the case for 
any physical operator of interest. The operators we are 
interested in can be classified as operators acting on the 
left dot or operators acting on the right dotcS, namely: 



(Ol)qi, = 1,Qh = 1 



r^c(l,0) 



'^c(l,l) 



{Oi 



\0,1 



^c(0,0) 



^c(l,l) 



(A9) 



We can conclude from this analysis that physical op- 
erators on the left and right sides have to be of or- 
der 0{e-'P^^) + 0(e-*^^^e-*''^«) and 0{e~'l^^'') + 



0{e 



), respectively. From now on we denote 



the order of the operators as 0(1, 0) -I- 0(1, 1) (left oper- 
ators) and 0(0, 1) -I- 0(1, 1) (right operators). 
Eqs. (A9) can be rewritten as: 



{Ol)ql=i,Qr=i 

(Ofl)Qi = l,Qjj = l 



^c(l,0) 
^c(l,l) 



\0,1 



^c(m,A .1.11 

^c(0,0) 



^c(l,0) 



{On) 



(AlO) 



The coefficients and can be identified with 

the right and left normalization factors in the absence of 

inter-dot honnintr i e ^^"^"'"^ = J- and = J- 

iniu aoi uoppmg, i.c, ZcCd.l) — Zn Zc{\.ff) ~ Zl 

and can be obtained from the left and right canonical 

partition functions of two independent single impurity 



problems at different chemical potentials /i^ and /^i^, re- 
spectively (see Refs. |lj,^,^. This way, the physical 
operators in the constrained ensemble can be written as: 



(Ol)qi,=i,q„=i 



^c(l,0) 



\0,1 



Zl 



(All) 



Eq. (All) is the central result of this section. Left 
(right) physical operators in the restricted ensemble 
with Ql = 1, Qb, — 1 contain two terms: 

i) the coefficient of the term of order 0(e~*'''^^('^' ) 
in the operator evaluated in the grand canonical ensem- 
ble plus 



a) the coefficient of the term of order 0(e 



in the operator evaluated in the grand canonical ensem- 
ble divided by the normalization factor of the right (left) 
problem without inter-dot coupling. 



The new normalization factors. 



Zc{0,l) 



and 



Zc(l,0) 



Zc(l,l) " " Zc(l,l) 

be obtained from the identities {Ql)ql=i,Qr=i = 1 

(<3-R>Qi.=1,Qh=1 = 1- 

Now, we apply the previous projection procedure to 
the selfenergies of Eqs. (10) and (11). The projection of 
the selfenergies can be accomplished in three basic steps. 
In a first step, we follow Langreth and Nordlander in 
Ref. ^ (see also Ref. |l^): since the Dyson equations for 
^L{R) ^^'^ ^L(R} *^ontain either G^^^^ or B^^^^ in every 
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term, the selfenergies that muhiply these quantities must projected out. As a result we obtain from Eqs. (10) and 
have all terms proportional to G^,„s or B^,„, or higher (11) the following selfenergies {t = t — t'): 



ni(«)(r) 



Similarly, 




+ *(^)'[GJi(L),.(r) - G?e(^),.(r)]B< ^)(-r)| X 
[^T.^Lina--)+^i'frT.^B^HiL^r)~B^^^^^^^ (A12) 

a a ) 



V 



, . L(R) 



FIG. 17: Self-energy contribution of order 0{e~^^^) to the 
conduction electron propagator which is projected out by the 
constraint. The inclusion of this selfenergy contribution to the 
conduction electron propagator would give unwanted contri- 
butions of order 0(1, 0) (left side) and 0(0, 1) (right side). 



Where we have emphasized in our notation the structure 
i^Kernell + Kernel2^ x propagator. In kernel 1, the 

conduction electron propagators K are defined in terms 
of the Fourier transforms of the bare conduction elec- 
tron propagators (namely, without dot-lead coupling) as: 

K|<.(e) = 2^Efc„^'^(e-efcJ/l(e), where /<(e) = 
is the Fermi function and /> (e) = 1 - /< (e) 
(see Ref. ^ and below). This way, the kernel is 0(0, 0). 
The Green's functions within the kernel2 part, namely 
G,B, do not include inter-dot hopping meaning that the 
kernel is 0(0, 1) for the left part and 0(1, 0) for the right 
one. This previous projection in the kernels is completely 



equivalent to the projection one does in the single impu- 
rity problem: with the same kind of arguments one ne- 
glects terms of order 0{e~'^^^) in the conduction electron 
propagator which, in principle, is a full propagator to be 
calculated in the presence of slave fermions and bosons. 
Self-energy corrections to the lead electron propagators 
like the bubble diagram in Fig. 17 are thrown away in 
the single impurity case (and also, of course, in our case). 
As a consequence of this projection, one always works 
with bare conduction electron propagators, which, again, 
is not what one gets initially from the unprojected NCA 
equations. 

In a second step we project out unwanted contribu- 
tions from the propagators multiplying the kernels. As 
we mentioned previously, the left (right) kernel is of or- 
der 0(0,0) -f 0(0,1) (0(0,0) + 0(1,0)), meaning that 
the retarded propagators multiplying these kernels Eq. 
(A12) should be of order 0(0,0), namely bare prop- 
agators with respect to interdot. The corresponding 
lesser propagators in Eq. (A13) therefore contribute with 
0(1, 0) (left operators) and 0(0, 1) (right operators) giv- 
ing the correct order: 0(1,0) + 0(1,1) for the left op- 
erators and 0(0, 1) -I- 0(1, 1) for the right ones. Finally, 
according to Eq. (All) the 0(1, 1) contributions should 
be normalized by Zj^ and respectively. The final set 
of projected selfenergies is then: 



ni(H)(r) = 



\ (7 (7 



X B 



L{R) 



G 



R(L),a 



i-r)' 



^R(L) 



^Gl^R)Ar). (A14) 



18 



Which correspond to Eqs. (19) and (20) used in the mam text. 
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Here, we focus on a projection procedure which is valid for 
diagonal (in the left/right index sense) operators. As we 
mentioned in section IVA, we can construct a fully self- 



20 

consistent NCA theory to order 0{1/N) by just including 
diagonal propagators/self-energies. 



